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ABSTRACT 

Context. The convective envelopes of cool main-sequence stars harbour magnetic fields with a complex global and local structure. 
These fields affect the near-surface convection and the outer stellar atmospheres in many ways and are responsible for the observable 
magnetic activity of stars. 

Aims. Our aim is to understand the local structure in unipolar regions with moderate average magnetic flux density. These correspond 
to plage regions covering a substantial fraction of the surface of the Sun (and likely also the surface of other Sun-like stars) during 
periods of high magnetic activity. 

Methods. We analyse the results of 18 local-box magnetohydrodynamics simulations covering the upper layers of the convection 
zones and the photospheres of cool main-sequence stars of spectral types F to early M. The average vertical field in these simulations 
ranges from 20 to 500 G. 

Results. We find a substantial variation of the properties of the surface magnetoconvection between main-sequence stars of different 
spectral types. As a consequence of a reduced efficiency of the convective collapse of flux tubes, M dwarfs lack bright magnetic 
structures in unipolar regions of moderate field strength. The spatial correlation between velocity and the magnetic field as well as the 
lifetime of magnetic structures and their sizes relative to the granules vary significantly along the model sequence of stellar types. 

Key words, magnetic fields - magnetohydrodynamics (MHD) - methods: numerical - stars: activity - stars: atmospheres - stars: 
low-mass 


1. Introduction 

Cool main-sequence stars have outer convection zones (spec¬ 
tral types F to about M4) or are fully convective (spectral types 
later than M4). The turbulent flows in the electrically conduct¬ 
ing plasma along with the stellar (differential) rotation support a 
hyd romagnetic dynamo that constantly generates magnetic flux 
(see ICharbonneaul[20 1 OL for a review). A fraction of this mag¬ 
netic flux becomes buoyant and reaches the near-surface layers 
where it is c ontinuously restructured by the convect ive flows (for 
reviews, see lFanll2009l:ISteinll2Q 1 2l: ISchussleill2Q 1 3l) . At the same 
time, the magnetic field strongly affects the stellar near-surface 
layers. As the solar example illustrates , the magnetic field causes 
a variety of activity phenomena (see ISchriiver & Zwaanll20dM 
for a review). While many of these (e. g. flares and coronal mass 
ejections) are located in the higher atmospheric layers and in the 
interplanetary space, the solar magnetic field also has an impact 
on the photosphere and the convection zone. Sunspots, pores, 
and faculae are well-studied photospheric phenomena that are 
caused by strong (1.5-3 kG) localised magnetic fields near the 
optical surface. Here, the magnetic field modifies the radiative 
and convective prop erties of the near-surface layers (see e. g. 
iRempel et al.ll20093) . 

Like the Sun, many stars exhibit magnetic activity in the 
for m of flares. X-ray emission , chromospheric emission lines, 
etc. (ISchriiver & Zwaanll2008l) . Lightcurve variability owing to 


spatial brightness variations (e. g. stars pots) is nowadays de - 
tected in essentially every cool star (e. g. iReinhold et al.l 120131) . 
Measurements of activity indicators (such as X-ray luminosity or 
the S-index) reveal that rapidly rotatin g (young) stars are m uch 
more magnetically active than the Sun (iReiners et al.ll20l4 and 
references therein). This also affects t he detection of pl ane- 
tary companions around these stars (e. g. I Jeffers et al.ll2013b and 
probably als o reduces the chanc e for habitability of the hosted 
planets (e. g. lVidotto et al.ll201^ . 


For stars other than the Sun, current observations cannot di¬ 
rectly spatially resolve activity-related phenomena in the pho¬ 
tospheres. While most activity studies rely on activity indica¬ 
tors, for some stars it is possible to measure the surface mag¬ 
netic field_with_spectroscopic and spectropolarimetric methods 
(see lRein^l20 1 21 for a review). For cool main-sequence stars, 
these measurements yield global magnetic held stre ngths rang¬ 
ing fr om a few Gauss ( from spectropolarimetr y, e. g. iPetit et all 
l2008l) up to > 3kG (iReiners & Basril [200% . However, it is 
known from the Sun that the magnetic held at the surface is 
highly structured and inhomogeneous on length scales down to 
the current resolution limit of solar observations (~ 100 km) and 
probably below. This local structure inevitably has a sig nificant 
effect on spectroscopic and spectropolarimetric data (IStenflol 
l2013l) . In order to better understand the global magnetic held 
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properties, knowledge of the local structure of the magnetised 
photosphere is therefore essential. 

Since the 1990s, comprehensive 3D radiative magne¬ 
tohydrodynamics (MHD) simulations have been used to 
simulate various magnetic phenomena in the solar near¬ 
surface layers 2 _ranging__^from_sim 2 le_jiet^wka^ plage re- 


dons (iNordlunc 

[& Stein! 119901; IVoderetal. l2005l) to Dores 

(ICameron et all 

2007b. sunsDots (iRemDel et all l2009ab. ac- 

five regions (ICheuns et al. l2010h. emergence of horizon- 

tal magnetic flux (IStein & Nordlundl 20061: IStein et al. 

|2^. 

and small-scale 

dynamo action (IVosler & Schiissler 

120071: 


iPietarila Graham et al.ll201Ql) . Until a few years ago, MHD sim¬ 
ulations of this kind were only available for the Sun. Recently, 
the first comprehensive 3D radiative MH D simulations for 
stars other than the Sun were presented ( Beeck et 'ZI 1201 U 
Wedemever et alJ l20 1 3t ISteineret~^l2013U20l4 iBeeck et all 
201 5W). 

In iBeeck et ^ d2013allR hereafter Paper I and Paper II), 
we discussed a sequence of non-magnetic (hydrodynamic) 3D 
simulations of near-surface convection with parameters close to 
main-sequence stars of spectral types ranging from F3 to M2. In 
Paper I, we analysed the overall structure of the convection and 
horizontally averaged stratifications. In Paper II, the granulation 
pattern in computed intensity images was studied and the effect 
of the 3D nature of the convection on limb darkening and on 
the profile of spectral lines was analysed. Here we continue this 
series of papers by considering the effects of a moderate mag¬ 
netic field on the near-surface convection in cool main-sequence 
stars. For each set of stellar parameters presented in Paper I and 
Paper II, we consider three MHD simulation runs with a unipolar 
average magnetic field of various strengths. A brief description 
of the simulation setup is outlined in Sect.O In Sect.O the distri¬ 
bution of the magnetic field near the optical surface is discussed. 
In Sects. |4] and [51 we study the structure of small-scale magnetic 
flux concentrations. In the following sections, the effect of the 
magnetic field on the heating of the upper photosphere (Sect. [6]) 
and on the convective flows (Sect. [ 7 ]) is analysed. A brief con¬ 
clusion follows in Sect. (H _ 

In a subsequent paper of the series (iBeeck et al.ll20I5al her- 
after Paper IV), we will consider the effects of the magnetic field 
on the limb darkening and on spectral lines. 


2. Simulation setup 

The simulations discussed in this pa per have been obtained with 
the 3p radiative MHD code MURaM (IVogleill2003l: IVogler et"^ 
l2QQ5l) . The version of the code used for the si mulations con¬ 
sidered here includes the modifications made by iRempel et al.l 
(l2009bh to enhance the efficiency of the magnetic runs. We limit 
the Alfven speed to Cmax = 60kms“^ to prevent exceedingly 
small time steps. This only affects the nearly force-free regions 
near the tops of the simulations boxes and does not considerably 
influence the results discussed in this paper and in Paper IV. The 
numeric al scheme is described in the appendix of iRempel et al.l 
(I2009bh . 

The radiative transfer module of the code utilises the opac¬ 
ity binning method (see discussion in Paper I). For the non¬ 
magnetic simulations considered here, opacity distrib ution func¬ 
tions provided by the ATLAS9 code (lKurucziri993l) were rear¬ 
ranged in four groups applying the r-sorting method (bin limits 
at logTref = 0, -2, -4; see discussion in Paper I). The magnetic 
simulations used the same binned opacities as the corresponding 
non-magnetic simulations. 


Six sets of stellar parameters (log g and Teff) roughly match¬ 
ing main-sequence stars of spectral types F3V, G2V, KOV, K5V, 
MOV, and M2V were considered. The horizontal and vertical 
sizes of the computational domains were scaled with the ex¬ 
pected granule sizes and pressure scale heights, respectively, to 
obtain similar setups for all sets of stellar parameters (for details, 
see Paper I). Snapshots of the non-magnetic simulations pre¬ 
sented in Paper I and Paper II were used as initial configuration 
for the mass density, the internal energy density, and the velocity 
field of the magnetic simulations. The initial magnetic field was 
unipolar, vertical with a uniform field strength of ^o- For each 
non-magnetic simulation, three simulation runs with = 20 G, 
100 G, and 500 G, respectively, and otherwise identical initial 
conditions were performed. In a transient restructuring phase, 
the flows advect the field into the converging downflow regions. 
After some time, the properties of magnetic field and flows be¬ 
come statistically stationary. The duration of the restructuring 
phase is of the order of the local convective turnover time, i. e. 
several minutes at the optical surface for all stellar parameters of 
our simulation sequence (and somewhat longer in the deeper lay¬ 
ers). The simulations were run considerably longer to make sure 
that a statistically stationary state was reached. Tabled] gives the 
runtime of the simulation ti (in simulated stellar time) at which 
the first snapshots were considered and the total runtime, 
of the simulations. The effective temperature, Teff, of the simu¬ 
lations is regulated only by the entropy density of the inflowing 
plasma at the bottom of the computational domain. This quantity 
was fixed and identical in all four runs for each set of stellar pa¬ 
rameters. Owing to the different (hot and cool) magnetic small- 
scale structures, there are slight differences in Teff between the 
different runs of the same star; values are given in Tabled) Values 
of Teff of the non-magnetic runs are almost identical to the ones 
of the 20 G runs and are given in Table 1 of Paper I. Overall the 
difference in Teff between runs of different is small. The rea¬ 
sons for the mostly positive deviations of the lOOG and 500 G 
runs from the non-magnetic and 20 G runs will be discussed in 
Sect.[6j 

The combination of three different values for with the six 
different sets of stellar parameters results in 18 magnetic simu¬ 
lations Q The six non-magnetic simulations presented in Paper I 
and Paper II serve as a reference, completing the grid of 24 sim¬ 
ulation runs. 


3. Magnetic field distribution at the optical surface 

Figure [T] shows the vertically emerging bolometric intensity of 
all three magnetic runs for the four spectral types F3V, G2V, 
KOV, and M2V. An analogous figure for the non-magnetic sim¬ 
ulations is given in Fig. 2 of Paper I. In all simulations, regions 
of either strongly enhanced or significantly reduced intensity ap¬ 
pear, mostly situated in the intergranular lanes. These regions co¬ 
incide with the locations of strong magnetic flux concentrations 
(see Fig. [2]). The relative number and size of dark and bright re¬ 
gions as well as their intensity contrast strongly depend on the 
spectral type and the amount of magnetic flux in the simulation 
box (i. e. on ^o)- This is discussed in some more detail in Sect.[4j 
The most striking difference appears between the M-star simu¬ 
lations on the one hand and all other simulations on the other: in 


^ Movies of the vertical velocity and the magnetic field strength at the 
optical surface as well as the bolometric intensity of the 18 magnetic 
simulations are provided as online material. They show the temporal 
evolution of the granulation during the time interval considered in this 
paper and Paper IV. 
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Fig. 1. Maps of the vertical bolometric intensity for twelve of the 18 magnetic simulation runs. For improved image contrast, the 
grey scale saturates at the values indicated by the colour scales on the right of each row. 


the F-, G- and K-stars simulations, only for Bq = 500 G a few 
dark features appear among the mostly bright magnetic struc¬ 
tures; in contrast, in the M2V star (and in the MOV star; not 
shown here), dark magnetic structures frequently form already 
in the 100 G runs and the intensity maps of the 500 G runs are 
dominated by large dark structures while there are es sentially no 
bright magnetic structures left fcf. lBeeck et al.ll2QTl1) . 

Figure [2] shows maps of the vertical component of the 
magnetic field at z = 0, the mean height level of the optical 


surface for the same selection of simulation runs and at the same 
time steps as Fig. [T] The strong concentration of the magnetic 
fiux in some of the convective downfiows leads to local field 
strengths of up to several kG in all stars, even in the case of 
Bo = 20 G. The maximum of the field strength at the optical 
surface depends only weakly on Bq, especially in the cooler stars 
of the model sequence: for instance, in the M2V star, it is 3.5 kG 
in the 20 G run and 4.5 kG in the 500 G run (more values given 
in Table [2]). 
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Fig. 2. Maps of the vertical component, of the magnetic field at z = 0, the average geometrical depth of the optical surface (same 

simulation runs and at the same time step as in Fig.[T}. 


As the typical field strength in the magnetic fiux concentra¬ 
tions does not strongly vary with the area fraction of these 
structures increases with increasing ^o- This is also illustrated by 
Fig.[3l which gives the histograms of the magnetic field modulus, 
B(tr = 1), at the optical surface. In the 500 G and 100 G runs, 
these histograms have two peaks (using logarithmically equidis¬ 
tant bins), one at relatively low values between 10 and 100 G (in 
most simulations, this peak falls between 30 and 50 G) and one at 
a high field strength of about 2 kG. The latter peak represents the 


typical field strength at the optical surface in the magnetic fiux 
concentrations and is denoted as ^strong in what follows. Values 
for ^strong obtained from a Gaussian fit to the histogram peak 
range between 1.4 kG and 2.4 kG (see Table [2]). ^strong slightly 
increases with Bq and towards the cooler end of our sequence. In 
the 20 G runs there is no strong-field peak, but the distributions 
of B(tr = I) still show an extended tail towards high values of 
the field strength. 
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Fig. 3. Histograms of the modulus of the magnetic field strength, B, at the optical surface for all 18 magnetic simulation runs (the 
histogram bin size is 0.04 dex). 


Table 1. Stellar parameters and runtime (in simulated minutes) 
of the simulations. 


Run (SpT-5o") 

logg[cgs] 

TeffiK] 


^end 

F3V-20G 

4.301 

6885 + 6 

111 

127 

F3V-100G 

4.301 

6911+8 

106 

121 

F3V-500G 

4.301 

7003 + 5 

116 

133 

G2V-20G 

4.438 

5779 + 9 

115 

133 

G2V-100G 

4.438 

5802 + 8 

108 

126 

G2V-500G 

4.438 

5864 + 9 

108 

126 

K0V-20G 

4.609 

4858 + 2 

102 

117 

KOV-IOOG 

4.609 

4878 + 4 

78 

91 

K0V-500G 

4.609 

4901 + 2 

74 

87 

K5V-20G 

4.699 

4376 + 2 

72 

83 

K5V-100G 

4.699 

4383 + 3 

61 

73 

K5V-500G 

4.699 

4402 + 2 

60 

72 

M0V-20G 

4.826 

3907 + 1 

83 

96 

MOV-IOOG 

4.826 

3909 + 1 

58 

69 

M0V-500G 

4.826 

3906 + 1 

51 

60 

M2V-20G 

4.826 

3691 + 1 

107 

124 

M2V-100G 

4.826 

3692 + 1 

59 

68 

M2V-500G 

4.826 

3679 + 1 

48 

56 


^ SpT: spectral type (of corresponding non-magnetic run); Bq\ aver¬ 
age vertical magnetic field strength 
^ ti : first time step considered in the analysis 
^ tend^ current end time of simulation 


The weak dependence of ^strong on is consistent with the 
rms values of B, which are largely independent of stellar pa¬ 
rameters. This supports a simple two-component model of the 
surface field: if one assumes a fraction / of the stellar surface 
to have a field strength ^strong and the remaining fraction 1 - / 
to have a field strength ^weak ~ 0 and further assumes the field 
to be mostly vertical and unipolar, one obtains / = ^o/^strong, 


Table 2. Properties of the magnetic field at the tr = 1 surface 
(all values in G; errors give temporal l-cr scatter) 


SpT" 

<S,) 

max(S(TR = 1)) 

V<52> 

D b 

-O strong 

F3V 

20 

2040 +160 

150 + 4 

— 

F3V 

100 

4250 + 550 

445 + 4 

1418 + 67 

F3V 

500 

5180 + 510 

1160+10 

1886 + 59 

G2V 

20 

2310+110 

178 + 6 

— 

G2V 

100 

3250 + 410 

412 + 10 

1702 + 74 

G2V 

500 

3550+ 150 

1022 + 11 

1990 + 18 

KOV 

20 

2890 +160 

175 + 7 

— 

KOV 

100 

2930 + 130 

425 + 3 

1824 + 24 

KOV 

500 

3390 +140 

1031+6 

2056 + 11 

K5V 

20 

2640+170 

176 + 4 

— 

K5V 

100 

2980 + 80 

418 + 6 

1823+48 

K5V 

500 

3550+140 

1043 + 3 

2099 +13 

MOV 

20 

3400 + 240 

185 + 5 

— 

MOV 

100 

3900 + 200 

436 + 4 

1917 + 41 

MOV 

500 

4200 + 200 

1075 + 3 

2326 +15 

M2V 

20 

3500 + 580 

182 + 3 

— 

M2V 

100 

4090 + 370 

440 + 8 

1994 + 85 

M2V 

500 

4450 + 320 

1067 + 7 

2352 + 57 


^ spectral type of the simulation 
^ peaks of the histograms of Fig. [3] 


from which for the rms = (BoSstrong)'^^- 

For ^strong = 2 kG one obtains an rms magnetic field of 200, 447, 
and 1000 G for = 20, 100, and 500 G, respectively, which is 
very close to the values actually found in the simulations (see 
Table O. Owing to the strong height dependence of the mag¬ 
netic field (cf. Fig.|4l), as well as to the local impact of the mag¬ 
netic field on the atmosphere structure, a two-component model 
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appears nonetheless inadequate to compute profiles of magneti¬ 
cally sensitive spectral lines (cf. Paper IV). 

4. Structure of magnetic flux concentrations 

Figure |4] shows vertical cuts through one of the magnetic fiux 
concentrations for each of the 100 G runs. The light-blue dashed 
curves give iso-contours of the magnetic field strength of 500, 
1000, and 1500 G, which indicate the locations of strong mag¬ 
netic fiux concentrations. In the fiux concentrations, the mag¬ 
netic pressure pmag = becomes substantial. As the at¬ 

mospheres are roughly in hydrostatic equilibrium, the gas pres¬ 
sure and the (mass) density are locally reduced in regions of 
strong magnetic field. The underlying grey-scale image gives 
the density distribution and confirms that the fiux concentrations 
are all under-dense with respect to the horizontal (iso-z) average 
stratification. However, the degree of evacuation strongly varies 
throughout the model sequence. The hottest two stars (F3V and 
G2V) show the strongest evacuation in the upper (optically thin) 
part of the magnetic structures, with a density of the order of 
only 10% of its horizontal mean value (in the first pressure scale 
height above the optical surface). In the K-dwarf simulations, the 
density in the magnetic structures reaches still about 30-50% of 
the horizontal average density and in the M-dwarf simulations, 
where evacuation is more confined to the upper layers, there is 
almost no evacuation near the optical surface. This effect is prob¬ 
ably a result of the much lower superadiabaticity in the cooler 
stars (see Paper I), which makes the convective collapse less ef¬ 
ficient (cf. iRaiaguru et al.ll200^ . The reduced efficiency of the 
convective collapse results in reduced ratio between the mag¬ 
netic pressure and the gas pressure (see Sect. O. 

The solid yellow contours in Fig. |4| outline the tr = I sur¬ 
face (= optical surface). In all simulations, magnetic fiux con¬ 
centrations cause a local depression in the optical surface, which 
is referred to as Wilson depression in analogy to the same phe¬ 
nomenon in sunspots. The depression is caused by the partial 
evacuation in the near-surface layers, which results in a lower 
absorption coefficient, kq, (and thus longer mean free path length 
of the photons). 

The red dash-dotted curves outline iso-T contours of 0.9, 1.0 
and 1.1 times the effective temperature. In all cases, the 1.1 
iso-contour has a lower geometrical height within the magnetic 
fiux concentrations than outside: at the same geometrical depth, 
all fiux concentrations are cool around z = 0. However, the op¬ 
tical surface depressions are in most cases larger than the de¬ 
pressions of the iso-T contours: at the same optical depth, most 
fiux concentrations shown here are hotter than the surroundings 
(with the notable exception of the M-star simulations). The con¬ 
tours of 1.0 and 0.9 Teff are mostly above the optical surface and 
show a more complicated structure. 

The optically thin plasma in the photospheric layers of the 
fiux concentrations is heated by various mechanisms including 
adiabatic heating in downflows (subadiabatic temperature gradi¬ 
ent), radiative heating, as well as viscous and Ohmic dissipation 
(see Sect. [6]). In the layers directly above the optical surface, the 
absorption coefficient is still high enough for radiative heating to 
be important. This leads to heating near inclined and concavely 
curved parts of the optical surface (side walls of the depressions). 
As a consequence, magnetic fiux concentrations with sizes com¬ 
parable to the Wilson depression have higher temperatures than 
the non-magnetic surroundings near the optical surface. In fiux 
concentrations with a large diameter (compared to their depth), 
this heating is only effective near the side walls. This explains 
why in the F-, G-, and K-type simulations the larger magnetic 


fiux concentrations become dark in their centres while smaller 
ones are bright. In the M-star simulations, the depressions in 
the optical surface are all very shallow owing to the low de¬ 
gree of evacuation in the upper part of the fiux concentration. 
Consequently, even a moderate amount of magnetic fiux gener¬ 
ates a structure with a large radius compared to its Wilson de¬ 
pression. Therefore, many magnetic fiux concentrations on M 
stars are dark. 

The dotted yellow contours in Fig. |4] give the tr = 0.01 
surface. Many photospheric spectral lines (such as the ones dis¬ 
cussed in Paper II and Paper IV) are most sensitive to the con¬ 
ditions around this optical depth. The magnetic field strength in 
the fiux concentrations and their filling factor may differ by a 
significant amount between the optical surface and this higher 
iso-TR surface. Note that the depressions in the tr = 0.01 sur¬ 
face are generally shallower than the depressions of the optical 
surface. This can partly be attributed to the heating mechanisms, 
which lead to a positive temperature deviation (with respect to 
the horizontal average) at this height (and above) and thus to a 
higher opacity, partially compensating the diminishing effect of 
the reduced density on the absorption coefficient. 

Figure [5] shows the dependence of the depth Wilson depres¬ 
sion on the spectral type of the simulation and on ^o- In the left 
panel of the figure, the depth is given in terms of km. By far the 
deepest Wilson depressions are observed in the F3 V star (about 
950 km in the 500 G run). The depth of the depressions decreases 
strongly along our simulation sequence towards cooler stars as 
well as with decreasing ^o- The dependence on is mainly 
caused by the different size of the structures: in most cases, the 
depressions are deeper in their central part than near their bound¬ 
aries (where the field strength is somewhat lower and where the 
side-wall heating can result in a significant increase in tempera¬ 
ture and thus opacity). As the fiux concentrations become larger 
at higher values of their area becomes, on average, larger 
in relation to their perimeter. Therefore, these “boundary” ef¬ 
fects become less important at larger resulting in a greater 
average depth of the Wilson depressions. In the right panel of 
Fig.B the depth of the Wilson depressions is given in terms 
of pressure scale height. Although the pressure scale height de¬ 
creases by more than one order of magnitude (from ~ 500 km in 
F3V to ~ 35 km in M2V, see Paper I), the decreasing trend of 
the depth of the depressions along the model sequence still re¬ 
mains. The depth of the depressions ranges from 0.17 (M2V, 

Bq = 20G) to \.5Hp (F3V, = 500G). Here, a threshold of 

5(tr = I) = 750G was applied to identify magnetic fiux con¬ 
centrations. Using different threshold values does not influence 
the results qualitatively, but has a quantitative impact, in particu¬ 
lar for the 20 G and 100 G runs: with a threshold of 500 G rather 
than 750 G, the depth of the Wilson depressions decreases by 
roughly 20%, 15%, and 5% for the 20 G, 100 G and 500 G runs, 
respectively. With a threshold of lOOOG rather than 750 G, the 
depth of the Wilson depression increases by roughly 15%, 10%, 
and 5% for the 20 G, 100 G and 500 G runs, respectively. The 
small impact on the 500 G runs confirms that the dependence on 
^0 is mostly due to the boundary effects discussed above. 

5. Balance of magnetic, thermal and turbulent 
pressures 

Owing to the Wilson depressions, the optical surface within the 
magnetic fiux concentrations is below its average level z = 0. 
As the magnetic field lines fan out in the upper parts of the fiux 
concentrations, the field strength is higher at the optical surface 
(where it defines ^strong) than at z = 0 (cf. Fig. (4]). The magnetic 
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Fig. 4. Vertical cuts through typical magnetic flux concentrations in the six simulations with Bq = 100 G. The grey scale indicates 
the density gl{g)z relative to the horizontal mean density; the solid (dotted) yellow curve represents the tr = 1 (tr = 0.01) surface; 
the dashed blue contours indicate iso-5 contours of 500, 1000, and 1500 G; the red dash-dotted curves indicate the T /Teff = 0.9, 
1.0, and 1.1 surfaces, x and z coordinates are to scale for each panel but differ by more than one order of magnitude for the different 
stellar models. 
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held strength {B{z = 0))int in the magnetic flux concentrations 
varies somewhat more strongly along our sequence of spectral 
types than ^strong, but still much less than one would expec t 
if the convective collapse for thin flux tubes (cf. ISpruitIIl976h 
were equally efficient for different stellar parameters: this would 
yield (5 (z = fl))int/(8^) — 5mag,int(^ = 0) OC /7gas,ext(^ — A)? 
i. e. the magnetic pressure inside the flux tubes (referred to as 
“internal” in what follows) should roughly scale with the ther¬ 
mal pressure outside of the flux tubes (referred to as “external” 


in what follows), as the end products of an efficient convec¬ 
tive collapse are strongly evacuated magnetic structures which 
are roughly in pressure equilibrium with the external gas pres¬ 
sure. Table [3] lists horizontal averages (at z = 0) of the inter¬ 
nal magnetic pressure, the external gas pressure, and the exter¬ 
nal “turbulent pressure” (twice the kinetic energy density) for 
the 500 G runs. Magnetic flux concentrations were defined as 
the area where 5(z = 0) > 750 G. The internal magnetic pres¬ 
sure in the K- and particularly in the M-star simulations does not 
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Fig. 5. Wilson depressions. Left panel: Difference in average geometrical depth of the tr = 1 surface between magnetic = 

1) > 750 G) and non-magnetic area (.S(tr = 1) < 750 G) as a measure for the depth of the Wilson depressions. Right panel: 
Downward shift of the optical surface within the magnetic flux concentrations in terms of pressure scale height (measured outside 
the flux concentrations). 


Table 3. Magnetic pressure (in the flux concentrations) and ex¬ 
ternal turbulent and thermal gas pressure at z = 0 in the 500 G 
runs in dyn cm“^. 


Simulation 

Pmag 

pv^ 


Pgas 


F3V 

7.41 

X 

lO"* 

3.21 

X 

To^ 

7.65 

X 


G2V 

9.39 

X 

lO"* 

3.87 

X 

lO'* 

1.46 

X 

10^ 

KOV 

1.10 

X 

10^ 

3.31 

X 

lO'* 

2.28 

X 

10^ 

K5V 

1.23 

X 

10^ 

3.49 

X 

lO'* 

3.15 

X 

10^ 

MOV 

1.47 

X 

10^ 

2.76 

X 

lO'* 

5.11 

X 

10^ 

M2V 

1.51 

X 

10^ 

2.00 

X 

lO'* 

5.99 

X 

10^ 


reach equipartition with the external gas pressure, as it approxi¬ 
mately does for the Sun and the F3 V star. In our model sequence 
(which roughly follows the main sequence through the Teff-log g 
plane), the decreasing efficiency of the convective collapse to¬ 
wards cooler stars is roughly compensated by the increasing ex¬ 
ternal gas pressure so that .^strong becomes almost independent of 
the stellar type in our sequence. 

Figure [6| shows 2D histograms of the magnetic held modu¬ 
lus versus the inclination, 0, of the held with respect to the sur¬ 
face (+90° = vertical; 0° = horizontal) in a layer covering half 
a pressure scale height both below and above z = 0. The mag¬ 
netic held in the near-surface layers falls into two components 
(cf. Fig. [3]): in regions with high held strength of a few kG (i. e. 
~ ^strong) the held is mostly close to vertical while weak-fleld 
regions show almost isotropic held inclinations with a weak ten¬ 
dency towards positive inclinations (i. e. the held direction of 
the dominant polarity). There is a dependence on spectral type: 
in the F3V-star simulation, the inclinations in the strong-fleld 
regime are less concentrated towards the vertical (with inclina¬ 
tions down to about 45°), while the held inclination in the M2V 
star is much closer to vertical (always above 60°) in this regime. 
In the latter simulation, the separation between weak and strong 
held is also more pronounced. Along the model sequence there 
is a smooth transition from the F-star-like to the M-star-like dis¬ 
tribution as illustrated by the other two spectral types shown in 
Fig. [6l This is related to the different balance of magnetic and 
turbulent pressures (see Table [3]): in the F3V star, the ratio be¬ 
tween internal magnetic pressure and external turbulent pressure 
is about 2.5 whereas it is 30 in the M2V star. This renders the 
magnetic flux concentrations in the F3V star much more sus¬ 


ceptible to the surrounding plasma motions, leading to a larger 
spread in held strength and orientation. In contrast, the magnetic 
flux concentrations in the M-star simulations are quite stationary 
and less strongly influenced by the surrounding convection. 

6. Heating of the upper photosphere 

Figure [7] gives the difference in {T)j{t^), the temperature aver¬ 
aged over iso-TR surfaces, between 500 G runs and non-magnetic 
runs, shown separately for up- and downflows. The local modifi¬ 
cation of the temperature in the flux concentrations shows in the 
average temperature profiles {T)j{t^). Interestingly, the impact 
of the magnetic held on the temperature structure is qualitatively 
and quantitatively similar in up- and downflows, although the 
magnetic flux concentrations are mostly located in downflows: 
the temperature is somewhat increased above the optical surface 
and decreased below the optical surface with the exception of 
the M stars, where the impact is smaller in amplitude and more 
complicated. The reduced temperature directly below the surface 
and the enhanced temperature above the surface imply a reduced 
average temperature gradient in the near-surface layers. This is 
a consequence of the somewhat more efficient radiative energy 
transport across the optical surface. Owing to the local depres¬ 
sions in the optical surface, the area of the radiating surface is 
increased. As a consequence, the radiative flux and hence the 
effective temperature of the atmospheres is somewhat higher in 
the magnetic runs. As the M stars only have very shallow Wilson 
depressions, they do not show a decreased temperature gradient 
around log tr = 0 in the magnetic runs. 

In most cases, the middle and upper photosphere (logTR < 
-2) is warmer in the magnetic runs than in the non-magnetic 
runs. While the heating of the lower photosphere above mag¬ 
netic flux concentrations is mainly due to radiation, the enhanced 
heating of the upper photosphere is caused by enhanced vis¬ 
cous dissipation and Ohmic heating in regions of substantial 
magnetic held. In the upper photosphere, the magnetic pressure 
dominates the thermal pressure. The magnetic field couples the 
plasma at this height to the magnetic flux concentrations below, 
which are subject to the interaction with the convective flows. 
This excites magneto-acoustic waves as well as torsional mo¬ 
tions, which are then dissipated at greater heights^For^larsim- 
ulat ions this effect has be en demonstrated by iMoll et aP (l20I2h 
and lShelvag et al.l (l2QI3l) . Observational evidence for magneto- 
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Fig. 6. 2D histograms of the modulus of the field strength versus the inclination of the field (0° = horizontal; +90° = vertical) for 
the 500 G runs of the F3V, G2V, KOV and M2V stars. 


acoustic waves in magnetic bri ght points in the Sun has been 
given bv iJafarzadeh et al.l (l2013l) . 

Figures [8] and [9] illustrate the heating of the upper photo¬ 
sphere. The figures represent horizontal cuts through the G2V-, 
KOV- and M2V-star simulations at a height of 4.6 pressure scale 
heights above the mean level of the optical surface. Figure [8] 
gives the non-magnetic runs of these stars, while Fig. [9] shows 
the lOOG runs. The upper row of plots illustrates the vertical 
velocity at that height. In the non-magnetic case (Fig. [8]), the 
velocity amplitude is relatively low and the horizontal scale of 
the structures is comparable to the granule size. In the lOOG 
runs, small-scale structures with high vertical velocities appear, 
which are associated with the magnetic structures (bottom row of 
Fig.O. The second row of the figures shows the modulus of the 
horizontal velocity at the same height.In the non-magnetic solar 
case, we find regions of high horizontal velocity that are usually 
accompanied by a shock front at an interface between up- and 
downfiows. In the cooler stars, the horizontal velocities are rela¬ 
tively low and no shocks appear. In the 100 G run, the differences 
between spectral types almost dissappear as all simulations show 
high-velocity horizontal flows with small horizontal length scale 
at the same places where the small-scale vertical velocity struc¬ 
tures appear. The third row of the figures shows the tempera¬ 
ture at the same height. In the non-magnetic solar simulation, 
the shocks are visible in the temperature as the localised viscous 
heating produces a high temperature at the shock fronts. In the 


non-magnetic runs of the cooler stars, the temperature fluctuates 
less at this level. In the magnetic simulations, the small-scale 
high-velocity structures correspond to regions of high tempera¬ 
ture fluctuations (and in most cases regions of higher tempera¬ 
ture). There are some localised maxima where the temperature 
is more than 1000 K (in the solar case even more than 2500 K) 
above the horizontal mean value. In Fig.[9l the fourth row shows 
the magnetic field strength at this height. The temperature max¬ 
ima are all located in regions with high magnetic field strength 
and often a significant gradient of the magnetic field. Although 
the general structure of the magnetic field is much smoother on 
the larger scales in this height, there is significant fine structure 
(e. g. current sheets and vortices), which is responsible for the 
strong heating on small scales. 

7. Impact on convective flows 

As described in Sect. O the velocity field is responsible for the 
very inhomogeneous structure of the magnetic field in the first 
place: horizontal outflows from the granules keep the magnetic 
flux in the downflow regions, while the up flows become nearly 
field-free. Owing to this process (in combination with the con¬ 
vective collapse), the field becomes locally very strong (several 
kG) and, consequently, the magnetic pressure can locally sur¬ 
pass the turbulent and thermal pressures (see Sect. IS. As a con¬ 
sequence, the magnetic field reacts back on the velocity field. 
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Fig. 7. Average temperature difference between the 500 G and non-magnetic runs in up- and downflows for the F3V, G2V, KOV, and 
M2V simulations. 


Figure [T0| displays maps of the magnetic held in a horizon¬ 
tal cut located 4.6 pressure scale heights below the optical sur¬ 
face for the 20 G and 500 G runs of the F3V, G2V, KOV and 
M2V stars. There are patches of both polarities. The dominat¬ 
ing polarity (here shown in red) is concentrated in the down¬ 
flows. Surrounding these are some regions of the opposite po¬ 
larity, which are produced by held lines caught in the overturn¬ 
ing motion at the transition between up- and downflows. In the 
right column of the figure, the vertical velocity at the same depth 
is shown for the 500 G run. As the mass supply in the mag¬ 
netic flux concentrations from horizontal inflows is suppressed 
by the strong held, the flux concentrations prevent the forma¬ 
tion of strong downdrafts below vertices of granules, which form 
in non-magnetic runs (see Paper I), and forces the downflows 
into a more network-like configuration around the magnetic flux 
concentrations and belo w unmagnetised intergranular lanes (cf. 
iNordlund & Steiniri99Ql) . 

Figure [TT] illustrates the effect of the magnetic held on the 
up- and downflow speeds averaged over surfaces of constant ge¬ 
ometrical depth z. As the magnetic held is concentrated in down¬ 
flows, the effect of the magnetic held is stronger here. Below the 
surface, the downflows are slower in the magnetic runs than in 
the non-magnetic runs. This effect becomes stronger for increas¬ 
ing values of ^o- In the higher photosphere, the magnetic runs 
often show a higher average downflow speed. This can be inter¬ 
preted as resulting from “convective collapse” events. The effect 
on the upflow speeds is qualitatively similar but smaller. It turns 
out that an average only over the weakly magnetised parts of the 
box (i. e. excluding the magnetic flux concentrations), shows no 
significant impact on the upflow speed near the optical surface 


whereas the downflow speed is still reduced by roughly 20% in 
the 500 G runs. This shows that the effect on the upflows at the 
surface is only due to the reduction of the area which is avail¬ 
able for normal convection, while the downflows are slower in 
the magnetic runs even in non-magnetic parts of the box. This 
is possibly a consequence of the more network-like structure of 
the downflows, which is more strongly affected by dissipation of 
kinetic energy than the efficient downdrafts in the non-magnetic 
runs. 

We also find a change in the rms of the horizontal velocities, 
which are smaller below the optical surface by 35 to 40% in the 
500 G runs and up to twice as high in the upper photosphere in 
the lOOG runs compared to the non-magnetic runs. The higher 
velocities in the upper photosphere can be explained by the mag¬ 
netic coupling to lower layers (see Fig.O. 

Figure [T^ shows 2D histograms of vertical velocity versus 
magnetic field strength in a layer of a width of half a pressure 
scale height both below and above z = 0. The (kinematic) weak- 
field regime has a broad distribution of vertical velocities as the 
upflows and a part of the downflows are only weakly magnetised. 
At field strengths of up to a few 100 G, downflows are favoured. 
This is at the end of the kinematic regime, where the magnetic 
held is passively advected into downflow regions and is not yet 
strong enough to react back on the flows. These would include 
magnetic flux concentrations during their (early) collapse phase. 
At high held strengths of a few kG, the distribution is rather nar¬ 
row around zero velocity with a tendency towards downflows; 
however, the differences between spectral types are most obvi¬ 
ous here: while the Sun and especially the F3V star still have a 
considerable spread in their velocities in the strong-held regime. 
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Fig. 8. Horizontal cuts through the boxes of the G2V, KOV, and M2V simulations (from left to right) without magnetic field. All cuts 
are situated such that {p)z = 0.01 po, i. e. about 4.6 pressure scale heights above the optical surface; the geometrical depth of this 
level is given in each panel. From top to bottom: Vertical velocity, horizontal speed, li^horl •= and temperature, T are 

given; see also Fig.[9l 


magnetic fiux concentrations are essentially at rest in the K and 
M stars. This probably is related to the different lifetimes of the 
magnetic fiux concentrations: while in the hotter stars, the in¬ 
dividual magnetic fiux concentrations are short-lived and being 
constantly rearranged and fragmented by the surrounding flows, 
the magnetic structures in the K- and M-dwarf simulations are 
much longer-lived, so that most of them are in a stationary phase 
at any given time. This can be understood in terms of the balance 
between the turbulent pressure, the thermal gas pressure, and 
the magnetic pressure in the fiux concentrations (see Table [3]). 
Owing to the much larger ratio of (internal) magnetic pressure 


to (external) turbulent pressure in the magnetic fiux concentra¬ 
tions of the cooler atmospheres of our sequence, the disruptive 
effects of the flows on the magnetic flux concentrations is much 
weaker. 

As the only source of information on the magnetic field of 
stars are spectra (more precisely, profiles of spectral lines sensi¬ 
tive to the Zeeman effect), the relation between the magnetic and 
velocity fields becomes a problem: active regions on a star have 
an rms velocity differing from quiet regions resulting in a dif¬ 
ferent velocity broadening of the spectral lines originating from 
differently magnetised surface components. Moreover, in our F- 
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row shows the horizontal field at the same height. 
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Fig. 10. Left two columns: maps of the vertical magnetic field component of the 20 and 500 G runs of the F3V, G2V, KOV, and M2V 
simulations at a height corresponding to 4.6 pressure scale heights below the mean level of the optical surface. The geometrical 
depth z of this level is specified in each panel. Right column: vertical component of the velocity of the 500 G runs at the same depth. 


and G-star simulations, the magnetic field is mostly located in 
downfiows, which adds a non-zero Doppler shift to the magnet¬ 
ically broadened line component. However, these effects are not 
taken into account at all when stellar magnetic fields are spec¬ 
troscopically determined. A detailed analysis of the effect of the 
combined magnetic and velocity fields on a few spectral lines is 
given in Paper IV of this series. 


8. Conclusion 

Solar observations show that a substantial fraction of the ob¬ 
served magnetic flux is distributed in “plage” regions with a 
moderate average field strength of up to a few 100 G, which lo¬ 
cally results in configurations very close to the one discussed 
in this paper. In these regions, magnetic bright points and fac- 
ulae as well as (micro-)pores are observed, which are similar 
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Fig. 11. Horizontal averages of the upflow {upper panels) and downflow {lower panels) speeds in the G2V {left panels) and M2V 
{rightpanels) simulations. 


to the structures visible in our and othe r comprehensive 3D 
RMHD siniulations of the Sun (see e. g. ICarlsson et al.l 120041 
iKeller et ^|2004 ICheung et al.ll2QQ8l) . 

Here, we investigate how the appearance of moderately mag¬ 
netised unipolar regions changes for a sequence of cool main- 
sequence stars (including the Sun). The major flndings are: 

- Th ere is a lack of bri ght magnetic structures on M dwarfs 
fcf. iBeeck et al.ll2QTl1) . Owing to the more efficient convec¬ 
tion (higher density near the optical surface, see Paper I) in 
these stars, the superadiabatic peak is less pronounced. This 
reduces the efficiency of field amplification by a co nvective 
collapse of thin flux tubes fcf. lRajaguru et al.|[2002h . The re¬ 
sulting Wilson depressions are shallow so that heating by 
radiation from the side walls becomes rather inefficient. 

- The magnetic field strength in small-scale magnetic flux con¬ 
centrations, .^strong, is almost independent of spectral type 
and of the total magnetic flux available. Along the main- 
sequence the thermal pressure at the optical surface increases 
with decreasing effective temperature. This compensates the 
lower efficiency of the convective collapse. 

- For all spectral types considered, the upper photosphere is 
heated similarly by viscous and ohmic dissapation above 
small-scale magnetic field concentrations. The magnetic 
field dynamically couples the upper layers (where /^mag ^ 
Pgas) to the deeper layers (where /^mag ^ /^gas)- 

- Convective flows are strongly modified within the magnetic 
flux concentrations. While the short-lived, smaller magnetic 
structures exhibit a strong downflow (convective collapse), 
the larger, long-lived structures are almost at rest as the in¬ 
flows into the structures are blocked by the strong held. 


We have shown that the spatial correlations between velocity 
held, thermodynamic variables, and the magnetic held qualita¬ 
tively and quantitatively differ between stellar types. The impact 
of this on the profiles of spectral lines in disc-integrated spectra, 
which are used for magnetic held measurements, is investigated 
in Paper IV of this series. 
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